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3 ■ Abstract. The fluxon solutions of a boundary problem for the sine- 

' ' Gordon equation (SGE) are investigated numerically in dependence on 

lO , the boundary conditions. Interconnection between fluxon and constant so- 

'"^ ' lutions is analyzed. Numerical results are discussed in context of the long 

^_^ ' Josephson junction model. 
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§^: 

c/3 , 1 Introduction 

-(-i ■ 

a The sine-Gordon equation is a nonlinear hyperbolic partial differential equation 
involving the d'Alembert operator and the sine of the unknovi^n function. The 

i-rt . equation reads 

Ph I utt - Ua^a; H- sin u = , (1) 

fj ' where u = u{x, t). It arises in differential geometry and various areas of physics, 

'— ^. including applications in relativistic field theory, Josephson junctions (JJs) or 

^^ ' mechanical transmission lines. The stack of coupled JJs describing by system of 

^ . coupled sine-Gorgon equations is investigating very intensively today |1I2I3J . 

^\ ' In the framework of the long JJs model, the dynamics of the magnetic flux 

0> . in a JJ of length 21 is described by the perturbed sine-Gordon equation: 

CN I (p^ja; -(/?«- a(pt = sin (^ - 7 , <>0, x&{-l,l) (2) 

t^ ■ 

^— V ■ with boundary conditions 

T-H \ ipx{±l,t) = he, (3) 

. . I where ip is the magnetic flux distribution, h^ - the external magnetic field, 7 - 

K^ ' the external current and a > - the dissipation coefficient. 

Jv>( I The aim of this work is a numerical investigation of the properties of the 

5_( ' static fluxon solutions of Eq.([2]) under the influence of the external magnetic 

Cu I field parameter /ig in ([3]) . Such solutions satisfy the following boundary problem 

- </3a:a; + sin(^-7 = 0, xe{-l;l), ipx{±l) = he. (4) 
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Here, we only consider the case of the JJ length 21 = 10 with zero external 
current 7 = 0. 



2 Numerical approach 

The static fluxon solutions of Eq.([5]) are obtained numerically, by solving of the 
boundary problem ^. The stability analysis is based on numerical solution of 
the following Sturm-Liouville problem [3] 

-i>xx + qix)ip = >^'ip, -ipxi±l)^0, q{x) = cos (p{x,p), p={l,he,-f). (5) 

In this approach, the minimal eigenvalue of Eq.([S]) Ao(p) > corresponds to a 
stable solution. In case Xq{p) < solution ip{x,p) is unstable. The case Xq{p) — 
indicates a bifurcation with respect to one of parameters p = {l,he,j). 

The numerical solving of Eq.Q is based of the continuous analog of Newton's 
method 0. At each Newtonian iteration the corresponding linearized problem 
is solved, on a uniform grids with 1025 number of nodes, using a three-point 
Numerov approximation of the fourth order accuracy [6] . 

For numerical solution of the Sturm-Liouville problem ^ we applied the 
standard three-point second order finite-difference formulae. First several eigen- 
values of the resulting algebraic three-diagonal eigenvalue problem are obtained 
by means of the standard EISPACK code. Details of numerical scheme are de- 
scribed in |7I8I9J for the double sine-Gordon equation. 

The known for h^ = solutions Mq and <P^ are numerically path-followed to 
non-zero positive he- At each zth step of the numerical continuation we analyze 
the stability of solution (p{x, he) and calculate the following physical character- 
istics: 

— full magnetic flux of the distribution Aip'^^^ip = ip{l, hi ) — '■p{~l, he); 

— quantity N denoted "number of fluxons" in [7] and determined as follows 

I 
N[^ix,hl:^)]^^Jcp{x,h<i^)dx. (6) 

Note, since each solution (f of Eq.(j3|) is defined with an accuracy 2fc7r (fc S Z) 
then the value N[ip] is also defined with accuracy 2k. The arbitrariness at the 
choice of integer number k can be used for the "concordance" of the value N 
with the value of the full magnetic fiux Atp according to the condition 

\N[ip]- Aip/27:\^mm. (7) 

Below, as in [71819111] . solutions ip with n — N where N satisfies Eq.© are 
denoted (/s". 

The crossing through the turning points in the numerical continuation (where 
the direction of the moving along the curve Aip[he) changes as we follow on the 
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new branch) was organized as in |10| . The turning points are identified with help 
of the relation that is tested at each ith step of numerical continuation: 



/..«-/i/~^' 



Z\(^(0 - Z\(y9(*^l) 



< e. 



(8) 



where e > is small known quantity. Note that ([5]) is a simple approximation of 
equality dh^/dAip = that is valid at the turning points of the curve Aip(he). In 
case we run into a turning point the sign of /ig-increment should be changed. At 
each step of the numerical continuation, the initial guess for Newtonian process 
was chosen in the form 



^(/^r^^) = ^{h^') + {hit'' hfi^) ■ '^^':(!!"r£^r'^ 



(9) 



that prevents the continuation from reversing to the previous branch of Aifijif,). 



3 Numerical results 
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Fig. 1. Internal magnetic field distri- Fig. 2. Internal magnetic field distri- 
bution ipx{x) associated with the state bution (px{x) associated with the state 
Mo for several values of the magnetic ^^ for several values of magnetic field 
field hp. hp. 



Two basic distributions are known at he = 0: the uniform Meissner solution 
Mo with N[Mo] = and the fluxon solution ^^ with N['^^] = 1 ([7]). Minimal 
eigenvalue Ag of Eq.(IS]) is negative for <P^ and positive for Mq. As we continue 
basic state Mq to he > Aq stays positive, i.e. the branch is stable until he — 2. 
In the <P^ case, the minimal eigenvalue Aq crosses zero at the point he = hi = 
0.054, i.e. the branch is unstable for < he < hi and stable for hi < he < 
2.098. Transformation of the internal magnetic field shape of basic solutions in 
dependence on he is shown on FigsHJ [51 
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Fig. 3. Dependence of the full magnetic flux A(p on the magnetic field he for 
fluxon distributions associated with Mq and <P^. Solid and dashed lines corre- 
spond, respectively, stable and unstable states. Light circles indicate the turning 
points, dark circles indicate the points where solution changes its stability. 
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Fig. 4. Dependence of the minimal 
eigenvalue Aq on the magnetic field 
he for the branch associated with 
<l>^. Light circles indicate the turning 
points, dark circles indicate the points 
of stability change. 



Fig. 5. Dependence of the minimal 
eigenvalue Aq on the magnetic field 
he for the branch associated with 
Mq. Light circles indicate the turning 
points, dark circles indicate the points 
of stability change. 



The ALp{he) branches associated with the basic solutions Mq and 'P^ are 
presented on Fig|31 It is seen, at some points (light circles in Fig|31) the curves 
A(p{he) turn back to another, upper, branches. When the Aip{he) curve turns to 
the left ("D"-point) the quantity N is increased to N + 2. So, the branch started 
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Fig. 6. Coexisting stable internal mag- Fig. 7. Coexisting stable internal mag- 
netic field distributions ipx at magnetic netic field distributions Lpx at magnetic 
field K ^ 1. field h, = 2. 



from the basic Mq solution a.t h^ = 0, joins the fluxons (stable and unstable) 
with the even N while the another branch (associated with the <P^ basic fluxon) 
connects fluxons (stable and unstable) with the odd N. 

The change of stability occurs at the points (marked by dark circles in FiglSj 
where the Ao(/ie) curve crosses zero, see Figs l4l5l The "d"- and "c" -turning 
points are indicated by the light circles. The "C" -turning points connect a pair 
of unstable solutions with the same number N: if" and ip". An increasing N to 
A^ -I- 2 is observed at the "D" -turning points (light circles). 

Thus, for < he < hi we have a single stable static distribution (associated 
with the basic solution Mq). For hi < h,, < her, her = 0.561 this distribution 
coexist with another one associated with the basic solution <P^ . An increasing of 
magnetic field he leads appearing, at her, (most left light circle in Fig [3]) a pair of 
(unstable) states ((p^, (fP'). As he is growing next, the stabilization of CfP' occurs 
(most left dark circle in Fig|31), i.e. for /le = 1 we have three stable distributions 
to be coexisting with unstable state (p^ , see FigjBl Further increasing he induces a 
creation, at each "c"-point, of additional pair (1^9", (^") with growing n, see Figl?) 
At the same time, the pairs ((^9" '^"+2 
disappearing at the "D" -points. 



i^"+2) with previous values n are sequentially 



4 Conclusions 



The detailed information on the variation of fluxon structure with external mag- 
netic field in long Josephson junction is very important for correct interpreta- 
tion of the experimental results. In this paper we investigated stationary fluxon 
solutions of Eq.(I5]),([2]) in dependence on the external magnetic fleld he- Our nu- 
merical technique allowed us to establish the interconnection between the basic 
solution Mq aX he ~ Q and the stationary distributions (ys" with even numbers 
n as well as the interconnection between basic state <?^ &i he ~ Q and (p" with 
odd numbers n. Coexistence of different stable n-fluxon distributions at different 
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values of external magnetic field he is been shown. We consider that predicted 
transformations of the stable fluxon distributions can be observed experimen- 
tally by investigation of the critical current in dependence of external magnetic 
field. 
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